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Abstract 

Inspired by classical work of Bel and Robinson, a natural purely algebraic con- 
struction of super-energy (s-e) tensors for arbitrary fields is presented, having good 
mathematical and physical properties. 

Remarkably, there appear quantities with mathematical characteristics of energy 
densities satisfying the dominant property, which provides s-e estimates useful for 
global results and helpful in other matters. 

For physical fields, higher order (super) n -energy tensors involving the field and 
its derivatives arise. In Special Relativity, they provide infinitely many conserved 
quantities. 

The interchange of s-e between different fields is shown. The discontinuity prop- 
agation law in Einstein-Maxwell fields is related to s-e tensors, providing quantities 
conserved along null hyper surf aces. Finally, conserved s-e currents are found for 
any minimally coupled scalar field whenever there is a Killing vector. 
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The importance of the Bel-Robinson and other super-energy (s-e) tensors JXL || 
is usually recognized, even though their physical interpretation remains somewhat 
obscure. Their mathematical usefulness is clear, a manifestation of which is their 
relevance in the proof of the global stability of Minkowski spacetime Q, or in the 
study of well-posed symmetric hyperbolic systems of differential equations including 
gravity (see, e.g. [|J]). In this paper, I show how to generalize these important 
mathematical properties to arbitrary fields and, more importantly, try to shed some 
light into the physical meaning of s-e by finding non-trivial conservation laws which 
involve two different fields. 

First, let us present the procedure to construct the s-e of any given tensor (5[. 
Consider any tensor i Ml ...^ m as an r-fold (ni, . . . ,n r )-form {n\ + ... + n r = m) by 
separating the m indices into r blocks, each containing ua (A = 1, . . . ,r) completely 
antisymmetric indices. This can always be done. Several examples are: F^ u = Fi^i 
is a simple (2)-form, while W pF^ u is a double (l,2)-form; the Riemann tensor R a px^ 
is a double symmetrical (2,2)-form (pairs can be interchanged); the Ricci tensor R^ u 
is a double symmetrical (l,l)-form, etcetera. Schematically, i^...^ will be denoted 
by i[m],...,[n r ] where [ua] indicates the A-ih. block with tia antisymmetrical indices. 
Duals can be defined by contracting each of the blocks with the volume element 
4-form, obtaining the tensors (obvious notation, signature -,+,+,+): 

t * ■ • • • j t * ^ i # * . . . . ^ t * *. 

[4-ni],...,[n r ] [ni],...,[4-n r ] [4-ni],[4-n 2 ],...,[n r ] [4-ni],...,[4-n r ] 

There are l+^^^+...+ ^^^=2 r tensors in this set (including i[ ni ] v .. 5 [ nr .]). 

Let us define the "semi-square" (£r ni i r^j x fr ni i r^i) by contracting all indices 
but one of each block in the product of t with itself (after reordering indices if 
necessary) 



(t x t) Xnil ,„x Hlr = l[ 




The s-e tensor of t is 2r-covariant and defined as half of the sum of the 2 r semi- 
squares constructed with tr ni i r^i and all its duals. Explicitly: 



2T Aim ... ArAlr {t} = [t [ni]r „ )[nr] x ^th],...,^]) + + 



t , . xi , * . (1) 

[i-ni],...,[i-n r ] [4-m],...,[4-»v]/ Xi f ii...X rf i r 

It can be proven that (jl]) is symmetric on each (Xa^a)-P^', if t[m],...,[n r ] 1S s y m_ 
metric in the change of two [n^-blocks, then (|lj) is symmetric with respect to the 
corresponding (Xa^a) -pairs; and (jl]) is traceless in any (Xa^a)-PSl^ with ha = 2. 
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For any timelike unit vector u, the s-e density oft relative to u is defined as 

W t (u) = T Xwi _.. Xr ^{t}u x ^ . . . u^u^ . 

Wt (u) is non-negative and satisfies 

{3u such that W t (u) = 0} <=> T AlMl ... ArMr {t} = <=> t^ m = 0. (2) 

Actually, Wt (u) is the sum of the squares |£ w ...^ m | 2 of all the components of t in 
any orthonormal basis {e^} with u = e*o. More importantly, the tensor ([!]) satisfies 
the dominant s-e property (DSEP) |], [7|, ||, ||, that is, for any future-pointing causal 
vectors ki , . . . , kir we have 

^Ai^i...A r ^t r {*}^l 1 ^2 1 • • • ^2r-1^2r — • 

This is equivalent to the "dominance" of To...o{t} in any orthonormal basis: 

W t (e ) = T ... {t} > |T AlMl ... ArMr {t}| , VAi,/xi, . . . , A r ,/i r . 

The proof of DSEP, which is one of the main properties of definition (|l]) , has been 
given by Bergqvist H in full generality using spinors. DSEP is a helpful tool to 
prove causal propagation of i Atl ,„ Atm (see e.g. ju], 0, and provides s-e estimates 
which may be useful in generalizing the theorems in || to the case with matter. 
Applying (Q) to the gravitational field, one constructs the tensors T a/3Xfl |-R[2][2]J 

and T a @ XfJ/ |C[2][2]| associated to the Riemann and Weyl curvatures, getting the 
classical Bel Q and Bel-Robinson (BR) tensors, respectively 

2 B a/3Xp = R apXa R l3 p p (7 + *R apXa *R l3 p p a +R* apXa R *> 3 p tl a +*R* apXa *R* f3 p t * (T , (3) 

£>a(3\ii _ £>(a(3)(\fi) _ gXfia/3 ^ Xfj, _ q-a(3Xfi _ rf(a/3\fl,) ^ <j~a Xfj, _ q 

These tensors have physical dimensions of L~ 4 , leading to two possible interpreta- 
tions: either as energy density "square" and references therein), or as energy 
density per unit surface. This second interpretation seems correct, as can be de- 
duced from several independent results [jl^, [l3| . Notably, the relationship 

E r = const. W r (u) \ r =o r 5 + 0(r 6 ) 

between the quasilocal energy E r of small 2-spheres of radius r in vacuum and the 
BR s-e density Wt {u) supports clearly this view (u orthogonal to the 2-sphere). 
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The Bel tensor satisfies 

VaB a{)\n = RP^JW + RP p ^J^P - ^g^R^J^P (4) 

where Jx P f3 = VxR^p — ^ P Rxp- Therefore, the Bel (as well as the BR) tensor is 
divergence- free in Einstein spaces (including vacuum): 

implying that T^^^grjxCu is divergence-free for any conformal Killing vectors £, ff, 
£. Other divergence-free tensors appear in Q], but they do not satisfy DSEP nor 
(|2j), and their lack of index symmetries prevents the existence of conserved currents. 

Let us now apply our general procedure to physical fields. Take any electromag- 
netic field F pu satisfying the source-free Maxwell equations: V p F u p = 0, Vu-F^i = 0. 
Using (0) one arrives at 

Tx,{F [2] } = Fx p Ef-±gx,F pcT F^ 

which is the standard divergence-free energy-momentum tensor, leading to con- 
served currents TP U {F[ 2 ]}Cu for any conformal Killing vector £. What about BR- like 
tensors? As a next step one can use the double (l,2)-form \7 a F pu as starting object 
to construct the corresponding tensor (|l|), which becomes 

T a( 3x^[i]F[2]} = V a F Xp V p Ff + V a F pp V p F x p - g a pV c ,Fx p V <T Ff - 

-\gx^ a F ap VpF^ + \g a pgx^rF ap V T F^ . (5) 

Chevreton's tensor Q is given simply by ^^{V^i^]} + Tx m p{V[i\F[ 2 ]}. The 
procedure can be continued building the (super) 2 -energy tensor T a px pT u associ- 
ated with the triple (l,l,2)-form V a V pFx p ., and so on, producing an infinite set of 
(super) n -energy tensors, one for each natural number n. The following fundamental 
result holds: the (super) 71 -energy (tensor) of the electromagnetic field vanishes at a 
point p iff the n th covariant derivative of F is zero at p. 

The divergences of (|5[) and higher order s-e tensors can be computed, and they 
vanish in flat spacetime providing infinitely many conserved quantities for the elec- 
tromagnetic field in Special Relativity. This result is not surprising, because if F pv 
satisfies Maxwell's equations in flat spacetime, then the 2-forms Vi ai i . . . \1 \ aa \F pv 
for fixed values of ax, ■ ■ ■ , a s so do, and their "energy-momentum" quantities give 
rise to the s-e tensors. In the full Einstein-Maxwell theory, there exists a divergence- 
free 4-index tensor || which, unfortunately, does not satisfy DSEP nor (|2|) and, due 
to its peculiar index symmetries, it does not provide conserved currents. 
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To explore the possibility of s-e exchange, let us consider the propagation of 
discontinuities |J. Assume there is an electromagnetic field propagating in a back- 
ground so that there is a (necessarily null) hypersurface of discontinuity E. Let 
us denote by [V]z the discontinuity of any object V across E. From the classical 
Hadamard theory one proves the existence of c p on E such that 

[•^/Lti/]j] — f^yfiv f^vC\ii TlpC^ — 

where n p is a null 1-form normal to E (notice that ft is in fact tangent to E because 
n^n^ = Ipll ), as well as the following propagation equation 

n^d^c 2 + c 2 (K + 2if>) = 0, c 2 = c M c M > 

where tj) is the factor appearing in n^V^n^ = ijjn 1 ' and K is the trace of the second 
fundamental form of E relative to n p (as E is null, such a second fundamental form 



is orthogonal to n and not extrinsic |16[| ). The above propagation equation ensures 
that if c p = at a spacelike 2-surface cut S of E, then c p = everywhere on E. For 
any conformal Killing vector £ it is then straightforward to prove that 



/ c 2 (n,cn 2 " 

Js 



is a conserved quantity along E, where u is the volume element 2-form of S. 
This can be related to the energy-momentum of the electromagnetic field because 
T[iv{F[2]} = c 2 n p n u if the background is empty. 

But what happens when [i ? At i/] s = 0? Then, T pu {F[ 2 ]} is continuous and s-e 
quantities must come in. Now there exist B pu and f p on E such that 

[Ra^x^ = 4n[ a B^n X ], B^ = B Ufl 2n^B p(3 + g^B^np = 0, 
[V a i^] s = 2n 0( n^/„], n p f p = 

and satisfying the propagation laws |l 

n^d^B 2 + B 2 (K + #) - 2n a F up BP T f T = 0, B 2 = B pv B» v > 0, (6) 
n»d p f 2 + f 2 {K + AiP) + 2n°F ap BP T f T = 0, f 2 = > . (7) 

Therefore, the discontinuities of the Riemann tensor and of V a F pv are a source for 



each other. In fact, from Q and (0) it is immediate that [17] 
^ t (B 2 + f 2 ) + (B 2 + f 2 ) (K + 4V) = 
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from where one finds the following conserved quantity along £ 

(8) 



J s 



which needs both the electromagnetic and gravitational contributions. The interest- 
ing point is that (||) is related to the tensors (|) and (||): if V a F pu vanishes at one 
side of S, then [B al3fl „}^ = 2B 2 n a n /3 n p n u and ^^{V^Fpj} = 2f 2 n a n f3 n p n u . 

To check whether or not this interplay of s-e quantities is generic, consider finally 
a scalar field (j) satisfying the Klein-Gordon equation V P V P 4> = m 2 (j) for mass m. 
Expression (|l]) for <f> would simply be T{(f)} = <ft 2 /2. However, the case of physical 
interest arises by constructing the tensor (||) associated to V^0: 

T X n{V [1]( t>} = V A </>V^ - ^g Xt y p 4>V p ^ 

which is the standard energy- momentum tensor of a massless scalar field. When 
m ytz 0, T{mcp} has the same physical dimensions than ?a m {V [i]</>} , so that one 
must combine them (using — g^ u which obviously satisfies DSEP) to obtain the 
whole energy-momentum tensor of the minimally coupled scalar field 

1 1 

Tx P = T Xfl {V ^cf)} + T{m(t>}(-gx p ) = Va^V^ - -g\ p V P ^ p 4> ~ ^rn 2 4> 2 g XlJ , ■ 

This procedure is systematic and the BR-like tensor for the scalar field is the sum of 
expression ([!]) for the double symmetric (l,l)-form V qS7 p(f> plus the terms coming 
from m. Expanding the duals one gets 

T al3 \ p = T aj3Xp {V y^V -r a/ 3{mV[i]0}5A M -r AiU {mV[i]0}5 a ,g +T{m 2 (j)}g a / 3 gx p = 
= 2V a V (A 0V M) V / # - g a p (VaV^V^V^ + to 2 V a 0V m </») - g Xp . (VaVW/jV^ 

+m 2 V^V^) + -g a pgx„ (v ct V p </>V ct V^ + 2m 2 V p </>V^ + m 4 2 ) (9) 

which satisfies DSEP, T a px p = T , (q/3)(A/j) = Tx pa p, and was previously found in 
Special Relativity [jl^]. Again one can construct (super) n -energy tensors associated 
to the higher derivatives \7 Pl . . . V Mn+1 (/>: the ( super ) n - energy (tensor) of the scalar 
field (m ^ 0) vanishes at a point p iff the covariant derivatives of <j) up to (n+l) th - 
order are zero atp. (If m = the result involves the (n+l)* ft -derivatives exclusively). 
Direct computation leads to 

V a T a ^ = 2V /3 V(a^ )p V^ - gxnWVpVpW^ 

-V^ (2VV (A ^ W + 5A^^ r VV>) (10) 

and similar longer formulae for (super) n -energy tensors, so that they are divergence- 
free in flat spacetime, providing infinitely many conserved quantities there. 
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The situation hitherto is that the Bel tensor (||) is divergence- free in vacuum, and 
the s-e tensor (||) of the scalar field is divergence-free in the absence of gravitation. 
The natural question arises: can these tensors be combined to produce a conserved 
quantity in General Relativity? To answer it, assume that the spacetime satisfies 
the Einstein-Klein-Gordon equations, so that 

iV = V^V.,0 + -m 2 2 ^ (11) 



and that £ is a Killing vector. Then £^V^0 = (if m ^ 0) |19| from where it also 
follows '0V M V 'u4> = 0. Using these results with @ and (0) one finds 

and from fll^) and (0) 

v Q (r^^CA^) = - (2VV A 0^ + la.V^) ^C A e 

so that finally 

Vai Q = o| r = (B a/3A " + . 



This also holds when m = for either case of £ being parallel or orthogonal to V0. 
The relevance of this result is that provides conserved s-e quantities (by means of the 
typical integration of j a and Gauss' theorem) proving the interchange of s-e between 
the gravitational and scalar fields, because neither B Q ^ ''^^ nor ^" £/?£a£m are 
divergence-free separately in general. These conservation laws make definition (|T|) 
of s-e tensors not only mathematically appealing but also physically very promising. 
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